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'  -  The  scheme  consists  of  two  major  computer  codes  which  in  turfi-consist  of 
a  large  number  of  subroutine  codes  and  of  empirical  relationships.  -First  of 
-all ,  >fhen  a  solar  flare  occurs,  six  flare  parameters  are  determined  as  the 
input  data  set  for  the  first  code  which  is  devised  to  show  the  simulated  prop¬ 
agation  of  solar  wind  disturbances  in  the  heliosphere  to  a  distance  of  2  AU. 
Thus,  one  can  determine  the  relative  location  of  the  propagating  disturbances 
with  the  earth's  position.  The  solar  wind  speed  -fV)'>-and  the  three  inter¬ 
planetary  magnetic  field  (IMF)  components  (Bx,  By,  Bz  or  B,  0,  <j>)"’are  then 
computed  as  a  function  of  time  at  the  earth's  location  or  any  other  desired 
(space  probe)  locations.  These  quantities  Cft-turn '  become  the  input  parameters 
for  the  second  major  code  which  computes  first^the  power  (e)-  of  the  solar 
wind-magnetosphere  dynamo  as  a  function  of  time.  The  power  thus  obtained  and 
the  three  IMF  components  can  be  used  to  compute  or  infer:  (1),'  the  predicted 
geometry  of  the  auroral  oval;  '-vti-i  )■*-/  the  cross-polar  cap  potential;  (444)-, > 

tTTe  two  geomagnetic  indices  A£  and  Dst^-  (4-v^/jthe  total  energy  injection  rate 
into  the  polar  ionosphere;  ,r-fv-}y  the  atmospheric  temperature,  etc.  r 

h\ 

At  the  present  time,  the  weakest  part  of  the  scheme  arises  from  the  fact 
that  causes  of  changes  of  the  IMF  Bz  (or  0)  component  are  not  well  under¬ 
stood.  Further,  several  empirical  relationships  are  incorporated  in  the 
present  scheme  in  the  second  major  code.  Thus,  there  are  still  a  number  of 
-solar  and  geomagnetic  storm  features  which  require  theoretical  studies  and 
proper  modeling  efforts  to  eventually  replace  the  empirical  relationships  by 
theo- etical  ones,  as  our  understanding  of  solar-terrestrial  physics 
a a vances . 
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1 .  Introduction 


1 . 1  Cause  of  Geomagnetic  Storms  and  Auroral  Activity 

Geomagnetic  storms  and  auroral  activity  are  caused  by  a  large-scale  elec¬ 
trical  discharge  process  surrounding  the  earth.  The  magnetic  fields  produced 
by  the  discharge  currents  are  identified  as  the  geomagnetic  storm  fields 
(namely,  a  magnetic  manifestation).  The  aurora  is  an  optical  manifestation  of 
the  discharge  process  which  is  caused  by  collisions  of  discharge-current 
carrying  electrons  with  the  polar  upper  atmospheric  atoms  and  molecules. 

1 . 2  Definition  of  the  Predict ion  of  Geomagnetic  Storms 

During  a  geomagnetic  storm,  an  intense  Van  Allen  belt  is  formed  near  the 
earth.  The  belt  carries  about  one  million  amperes  of  electric  currents,  and 
is  called  the  ring  current  belt  or  the  storm-time  Van  Allen  belt.  The  ring 
current  produces  a  southward  directed  field  in  low  latitudes  in  space  between 
the  earth's  surface  and  the  inner  boundary  of  the  belt.  The  magnetic  field  of 
the  ring  current  is  quantified  by  the  geomagnetic  index  Dst. 

A  pair  of  concentrated  currents,  called  the  auroral  electrojets,  are  also 
generated  in  the  polar  ionosphere.  The  total  current  intensity  of  the  elec¬ 
trojets  is  of  the  order  of  a  few  million  amperes.  They  cause  intense  and 
complicated  magnetic  fields  in  the  auroral  oval  and  its  vicinity.  The 
magnetic  fields  produced  by  this  current  system  are  quantified  by  the  geo¬ 
magnetic  index  AE.  The  AE  index  can  be  considered  as  a  complementary  index 
for  auroral  activity. 

The  geomagnetic  storm  intensity  is  thus  parameterized  and  quantified  by 
these  two  geomagnetic  indices,  AE  and  Dst.  Therefore,  a  modern  geomagnetic 
prediction  scheme  3hould  be  able  to  predict  AE  and  Dst  as  a  function  of  time  t 
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within  several  hours  of  a  major  solar  event;  note  that  other  complementary 
indices,  such  as  Kp  and  ap  indices,  are  derivable  from  the  AE  index. 

In  addition  to  the  geomagnetic  indices,  several  other  magnetoapher ic  and 
ionospheric  quantities  should  be  predicted  by  a  modern  storm  prediction 
scheme.  Some  of  them  are: 

(i)  Geometry  and  size  of  the  open  region  polar  cap  and  the  auroral 
oval . 

(ii)  Cross-polar  cap  potential  4>  _  . 

(iii)  Total  energy  injection  rate  into  the  polar  ionosphere, 
liv)  Increase  of  the  upper  atmospheric  temperature  AT^. 
iv)  Electron  density  profile  in  the  polar  ionosphere. 

1  •  >  i\sea  for  Numerical  Prediction  Schemes 

In  the  past,  the  geomagnetic  storm  prediction  has  entirely  depended  on 
statistical  results.  Two  examples  of  such  statistical  results  are  shown  in 
Figures  1.1  and  1.2  (Akasofu  and  Chapman,  1972).  Figure  1.1  shows  the  time 
lapse  from  flare  onset  to  storm  onset  (the  time  of  ssc).  One  can  see  that  the 
’.apse  time  is  about  *10  hours  on  the  average,  but  it  ranges  from  about  20  hours 
to  t> 0  hours.  Therefore,  there  is  a  great  uncertainty  in  predicting  the  onset 
time  of  a  geomagnetic  storm  after  a  flare.  Figure  1.2  shows  the  central 
meridian  dependence  of  the  maximum  Dst  value  attained  by  a  number  of  storms. 
One  can  see  that  there  is  a  clear  dependence  of  the  maximum  possible  Dst  on 
the  central  meridian  distance.  In  particular,  except  for  one  event,  all  major 
storms  [Dst  >  2U0Y  (=  nT)]  are  caused  by  flares  which  occurred  within  30°  of 
the  central  meridian  distance.  Such  a  tendency  is  important  in  predicting  the 
maximum  intensity  of  geomagnetic  storms.  However,  not  all  central  meridian 
flares  produce  large  Dst  storms.  In  Figure  1.2,  the  Dst  values  range  from  0 
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Histogram  of  the  transit  time  of  solar  flare- induced  disturbances  to  the  earth. 


to  H00  nT  within  the  central  meridian  range  of  ±20°,  so  that  one  can  provide 
only  the  largest  possible  Dst  value  for  a  given  flare  by  knowing  its  central 
meridian  distance.  For  this  reason,  there  is,  at  present,  no  way  to  predict 
the  time-development  of  the  Dst  index  f or  a  given  flare.  Therefore,  one  must 
develop  a  numerical  scheme  in  predicting  AE(t),  Dst(t),  the  geometry  of  the 
oval,  <tpc(t),  Uj  (t ) ,  ATm,  etc. 

1  . i)  Key  Element  in  the  First  Generation  Prediction  Scheme 

The  key  element  which  has  made  the  numerical  prediction  of  the  intensity 
of  geomagnetic  storms  and  auroral  activity  possible  is  that  we  have  succeeded 
in  identifying  the  solar  wind  quantities  which  are  related  to  the  power  t  for 
the  auroral  discharge  process  or  the  power  of  the  solar  wind-magnetosphere 
dynamo  (Perreault  and  Akasofu,  1978;  and  Akasofu,  1981).  It  is  given  by 

£  =  VB^sin4  (|)  l  2  (1  ) 

2  o 

where 

V  =  the  solar  wind  speed 

3  =  the  solar  wind  magnetic  field  magnitude 

9  =  the  polar  angle  of  the  solar  wind  magnetic  field  vector 

iQ  =  constant  (-7  earth  radii) 

A  tentative  criteria  for  most  systems  which  utilize  the  polar  ionosphere 
ar  e : 

d  1  g 

e  <  '0'  Mwatt  i < 1 0  erg/sec)  -  Normal  operation 

c  1 H 

c  -  bx10''  Mwatt  (-5  x  10  °  erg/sec)  -  Require  alert 
e  -  ’0^  Mwatt  (-10'^  erg/sec)  -  Prepare  for  system  failure 

e  >  Mwatt  9  x  10  ^  erg/sec)  -  System  failure 
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Tnerefore,  the  success  of  the  prediction  will  depend  on  tne  accuracy  of 
predicting  the  power  e(t)  as  a  f unction  of  time  within  several  hours  of  a 
major  solar  activity. 


2.  A_  First  Generation  Geomagnetic  Storm  Prediction  Scheme 

Figure  2.1  shows  our  first  generation  prediction  scheme  in  the  form  of  a 
block  diagram.  It  can  be  divided  into  two  parts,  A  and  B. 


Prediction  of  Disturbed  Solar  Wind  Quantities 


step  of  the  prediction  of  geomagnetic  storms  and  thus  of  the 
reduceu  to  the  prediction  of  solar  wind  quantities  V(t),  B(t) 
.t) ,  P..'t),  B  ( t ) ,  B . , '  t ) )  in  equation  (1).  Our  first  generation 

a  y  ^ 

I y  based  on  the  method  developed  by  Hakamada  and  Akasofu  (1982) 
...  M9 dp;.  This  method  provides  a  first  order  construction, 
sputi'iiiy ,  of  f iare-generated  shocks.  It  does  not  simulate  the 
thermodynamic  properties  that  can  only  be  found  from  the  MHD 
t none  qualities  are  not  essential  for  the  geomagnetic  storm 
r-  nave  already  been  a  number  of  simulation  studies  of  solar 
wmch  art,-  based  on  numerical  solutions  of  the  MHD  equations 
*  193i;  Gislason  et  al . ,  19 8^).  However,  a  full  three- 

2  simulation  .nas  not  been  developed  yet  for  our  purpose, 
i  is  ;o:sibie  to  ’calibrate’  our  kinematic  method  with  those 
two-dimensional  MHD  simulations  (Olmsted  and  Akasofu,  198b),  so 
should  be  reasonably  consistent  with  the  MHD  results.  In  the 


■ief  explanation  of  tne  codes  is  given, 


( i )  Basic  bac Kground  solar  wind 


This  code  sets  up  the  background  solar  wind  conditions  into  which 
disturbances  caused  by  a  major  sola,  event  will  be  introduced.  In  the 
present  scheme,  the  background  solar  wind  flow  can  be  specified  by  know¬ 
ing  the  location  of  the  magnetic  equator  in  heliographic  coordinates  and 
the  latitudinal  (magnetic)  variations  of  the  solar  wind  speed,  both  on 
the  so-called  'source  surface'  (a  spherical  surface  of  radius  of  rs  =  2.5 
solar  radii). 

Magnet i c  equator 

The  location  of  the  magnetic  equator  (or  the  so-called  neutral 
line)  on  the  source  surface  plays  a  crucial  role  in  modeling  the  back¬ 
ground  solar  wind  flow.  It  is  provided  by  the  Stanford  Solar  Observa¬ 
tory,  Stanford  University,  every  27  days,  namely  after  each  Carrington 
Rotation.  Therefore,  the  location  of  the  magnetic  equator  can  be  updated 
every  27~day  interval.  The  Stanford  data  also  provides  the  magnetic 
field  intensity  B  on  the  source  surface  and  the  intensity  Bp  at  the 
magnetic  pole. 

Latitudinal  variations  of  the  solar  wind  speed 

The  latitudinal  variation  of  the  solar  wind  can  be  approximately 
expressed  by: 


v  =  v0  +  (v,  -  vQ)  (1-cos0tXm) 
where  Am  =  sin-1  (B/B  ) 

The  three  parameters  a,  vQ  and  v-j  in  the  nth  year  of  the  sunspot 
cycle  are  tentatively  given  in  the  following  table  (Fry,  1985).  In  the 
scheme,  the  three  parameter  values  are  automatically  selected  by  speci¬ 
fying  the  year  in  A.D. 


r 

1st 

2nd 

3rd 

Uh 

5th 

6th 

7th 

8th 

9th 

10th 

11th 

32 

16 

8 

JJ 

3 

3 

JJ 

6 

12 

20 

32 

270 

270 

270 

270 

270 

270 

270 

270 

270 

270 

270 

700 

680 

670 

650 

650 

650 

660 

670 

680 

690 

700 

Background  code 

The  background  code  provides  the  solar  wind  speed  V,  density  n  and 
the  three  components  Bx,  B  and  Bz  (or  the  magnitude  B,  latitude  angle  0 
and  azimuth  angle  <p ;  9  =  90°  -  0  at  the  location  of  the  earth  for  a  27- 
day  period  and/or  at  any  other  location  in  the  heliosphere  for  space 
probes.  Two  displays  are  prepared  for  the  background  solar  win!  informa- 


OUTPUT  1 


The  first  is  the  location  of  the  magnetic  equator  in  Carrington 
longitude-latitude  coordinates,  together  with  some  iso-B  contours.  As  an 
example,  it  is  assumed  that  the  magnetic  equator  in  this  coordinate 
system  is  given  by  a  sinusoidal  curve  during  the  Carrington  Rotation  l8.-;6 
(November  23  -  December  19,  1990).  The  projected  location  of  the  earth 
on  the  source  surface  on  December  10,  1990  is  indicated  by  *. 

OUTPUT  2 

The  second  is  the  predicted  27-day  variation  of  the  solar  wind 
speed  and  the  density,  the  IMF  magnitude  B  and  the  two  angles  0  and  <*>  at 
the  earth's  location,  without  flares  or  other  solar  activities,  for  the 
same  Carrington  Rotation  period,  November  23  -  December  19,  1990.  Note 
that  since  the  magnetic  equator  is  sinusoidal  in  the  Carrington  coordin¬ 
ate  system  in  this  illustrative  example,  two  maxima  of  the  solar  wind 
speed  are  observed  at  the  earth.  Note  also  a  sharp  change  of  the  angle  <t> 
either  from  135°  to  315°  or  from  315°  to  135°,  indicating  the  crossing  of 
the  solar  current  sheet  (or  the  so-called  'sector-boundaries'). 

For  an  accurate  storm  prediction,  Output  1  and  Output  2  should  be 
continuously  updated,  since  flare-induced  disturbances  wiLl  be  superposed 
on  the  background  solar  wind  flow. 

( i i )  Sol a r  Flare  Code 

Solar  observatories  provide  the  following  information  on  a  solar 

flare: 


Start  time 

Time  of  maximum  flare  intensity 


• 

Latitude  (N,  S) 

• 

Central  meridian 

distance  (E,  W) 

• 

Importance  (0,  1 

,  2,  3,  H) 

• 

End  time 

In 

the  flare  code. 

each  flare  is  described  by  the  following  six 

parameters . 

• 

Start  time 

Tp:  Flare  start  time  minus  simulation  start 

time,  in  hours. 

. 

Latitude 

<j>p:  Flare  latitude  (N  =  +,  S  =  -). 

• 

Longitude 

Flare  longitude  (W  =  +,  E  =  1) 
earth’s  longitude. 

• 

Initial  speed 

Vp:  Flare  speed  at  source  surface  at  maximum 

epoch .  See  Table  1 . 

Area  <jp: 

Aerial  extent  of  flare,  standard  deviation  for 
gaussian  spatial  extent  on  source  surface. 

See  following  Table. 

• 

Duration 

ip:  Flare  time  constant. 

The  first  three 

parameters  Tp,  <f>p  and  ip  are  the  same  as  the 

observed  start  time, 

latitude  and  central  meridian  distance.  Based  on 

the  observed  Importance  unc  Brilliance,  Vp  ana  Op  are  chosen  from  the 
following  tentative  conversion  table.  In  the  scheme,  as  one  specifies 
the  flare  intensity,  this  procedure  is  automatically  accomplished.  The 
flare  time  constant  Tp  is  tentatively  taken  to  be  ten  times  of  flare 
duration  (see  Section  5(B)). 

Conversion  from  flare  Importance  and  Brilliance  to  VK  and  Op. 


Importance/Brilliance  _Vp  (km/sec)  _0p(°) 

OF  200  20 
ON  300  30 
0B  H00  H0 

IF  H00  30 
IN  600  HO 
1 8  800  50 

2F  500  30 
2N  800  50 
2B  1200  70 

3F  700  30 
3N  1100  50 
3B  1500  70 

HF  800  HO 
HN  1200  60 
HS  1600  80 


ijjr_t  B:  Predict  iori_of  Geomagnetic  l) i  st urbane es 

The  second  part  (B)  of  the  prediction  scheme  takes  the  outputs  of  the 
first  part  (A)  as  the  input.  As  the  solar  wind  speed  (V)  and  the  three  inter¬ 
planetary  magnetic  components  (Bx,  By,  Bz)  become  available  as  a  function  of 
time,  it  is  possible  to  estimate  the  auroral  discharge  power  e  as  a  function 
of  time  from  equation  (1).  On  the  basis  of  V,  _B(BX,  By,  Bz)  and  z  thus 
determined,  the  following  quantities  can  be  computed. 

(i)  Open  region  (auroral  oval)  geometry 

The  oval  geometry  eoae  takes  the  power  z  as  the  input  and  deter¬ 
mines  the  geometry  of  the  poleward  boundary  of  the  auroral  oval  in  mag¬ 
netic  latitude  -  local  time  coordinates.  This  code  is  based  on  a  study 
of  the  dependence  of  the  polar  cap  (the  open  field  line  region)  on  the 
IMF  by  Akasofu  ana  Roederer  (198H). 

(ii)  Gross-polar  potential 

The  cross-polar  cap  potential  <t>pC  i-  founa  to  be  proportional  to  /j7 
(Re  iff  et  al . ,  1981);  it  is  given  by  4>  =  30  +  0.006/F7  where  z'  = 

ii  p 

VB'sin  (6/2)  (10  nT  km/sec).  Combining  this  empirical  result  with  the 
auroral  oval  geometry  code,  one  can  also  infer  the  distribution  of  the 
electric  field  in  the  polar  regions.  We  also  have  the  option  of  using 
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the  empirical  relationship  between  <f>  and  AE,  determined  by  Ahn  et  al  . 
(1984);  it  is  given  by  4>  (KV)  =  36  +  0.089  AE(nT)  or  0  (K V)  =  ^.2  x 

”17  /  /  \  r 

10  '  e(erg/sec) . 

( i i i  )  AE,  Dst  prediction 

This  code  takes  the  power  e  as  the  input  and  predicts  the  AE  and 
Dst  indices.  The  AE  index  is  computed  on  the  basis  of  an  empirical 
relationship  between  e  and  AE,  AE(nT)  =  600[tanh (log1 qC- 1 8. 9 ) * 
650.  The  Dst(nT)  index  is  computed  from  equation  4  x  10^®  (8Dst/cit  + 
Dst/iR)  =  0.7e,  where  tr  is  the  lifetime  of  ring  current  particle.-', 
(Akasofu,  1982). 


( i v )  Magnetospheric  convect ion 

This  code  takes  the  power  e  as  the  input  and  determines  the 
convection  pattern  of  magnetospheric  plasma  over  the  entire  polar  iono¬ 
sphere.  It  is  assumed  that  the  cross-polar  potential  drop  4-  is  npnlieu 
to  the  open  field  line  region  determined  in  (ii). 


(v)  Total  energy  injection  rate  Uj  into  the  polar  ionosphere 


This  estimate  is  based  on  the  formula,  Uj  (watt)  =  2.9  x  10°  AE(nT) 
(or  2.9  x  lO1*  AE(nT)  erg/sec);  see  Ahn  et  al .  (  1984). 


(vi)  Increased  temperature  in  the  thermosphere 

This  estimate  is  based  on  the  formula,  AT_(K)  =  ex  -300  where  x  = 

17  00 

0.45  i,n  U(erg)/10  '  +  0.3)  +  5.0.  This  empirical  relation  w-s  obtained 
by  PrOlss  (1985) . 


3.  A  Cause  of  the  IMF  B_  Variations 

We  represent  schematically  the  sequence  of  magnetic  field  configurations 
which  might  result  from  a  particular  type  of  coronal  disturbance.  The  corona 
is  assumed  to  initially  have  a  dipolar  magnetic  field  configuration.  Figure 
3.1  shows  schematically  a  series  of  field  configurations  at  successive  time 
steps.  As  a  coronal  transient  occurs  in  the  lower  solar  atmosphere,  the 
magnetic  field  lines  are  forced  to  move  Doleward.  As  a  result,  the  inter¬ 
secting  points  of  the  field  lines  on  the  source  surface,  at  rs  =  2.5  solar 
radii,  move  poleward  This  induces  a  polar  component  in  the  field  line  orien¬ 
tation.  outward  from  the  source  surface,  and  a  polar  IMF  vector  component 


results.  The  distortion  of  the  field  lines  propagates  outward. 


As  the  dis- 


turbance  moves  outward  past  the  source  surface,  the  field  lines  gradually 
return  to  their  original  latitudes  on  the  source  surface.  Since  the  field 
lines  are  imbedded  in  the  outflowing  plasma,  the  distortion  continues  to  prop¬ 
agate  outward  beyond  the  source  surface. 

Figures  3.2  and  3.3  represent  the  meridional  plane  view  of  the  solar 
source  surface  in  the  case  of  a  poleward  motion  of  the  field  lines.  The 
kinematic  code  traces  the  field  line  as  follows.  At  time  t  =  t,,  the  sun 
emits  a  fluid  parcel,  or  "particle"  (position  1).  At  a  later  time,  t  =  t?, 
the  emitted  particle  has  moved  radially  outward  to  position  2.  The  source 
region  (which  emitted  the  particle)  has  moved  to  position  3  by  the  time  t  = 
t2-  The  region  emits  a  second  particle  at  position  3.  Particles  at  positions 
2  and  3  are  connected  by  the  same  field  line.  This  situation  is  similar  to 
the  streamline  of  water  formed  by  suddenly  moving  a  garden  hose  upward.  The 
resulting  kink  in  the  streamline  (and  the  magnetic  field  line)  propagates 
outward  at  the  fluid  flow  speed.  The  meridional  velocity  must  be  comparable 
to  the  radial  velocity  to  induce  a  0  angle  of  about  45°. 

_  ^ 

For  a  solar  wind  speed  on  the  source  surface  of  aoout  50  km  s  1  ,  the 

polar  motion  required  to  allow  Bz  comparable  to  Bx  in  the  IMF  is  also  about  50 
km  s'1.  This  is  because  even  though  the  radial  speed  vp  is  increasing  outward 

from  the  sun  due  to  acceleration  of  the  coronal  plasma,  vp  is  also  increasing 

outward  at  all  latitudes.  In  this  case,  the  meridional  field  line  kink  will 
experience  a  steepening  in  B  angle  as  it  propagates  outward,  since  the 
background  solar  wind  velocity  increases  away  from  the  magnetic  neutral 

line . 

The  field  line  location  on  the  source  surface  i3  reprosL-ut.-d  by  its 
magnetic  latitude,  which  is  now  allowed  to  vary  according  to  the  equation, 
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Fig.  3-3.  Top:  Schematic  showing  the  steepening  of  IMF  field  line  config¬ 
uration  as  the  kink  propagates  outward.  Velocity  V 2  >  V ^ .  Bottom:  Schematic 
snapshots  of  a  field  line  segment  at  successive  times  as  it  propagates  out- 
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X/mo  =  magnetic  latitude  for  the  dipole  case,  or  approximately,  X^ 

=  xsin (<p-4>o ) . 

The  background  magnetic  field  configuration,  Amo,  in  these  cases  is  that 
of  a  dipole  field  with  a  tilt  of  x  =  20°  with  respect  to  the  rotational 
axis.  At  observation  longitude  =  0°,  the  6  variation  is  similar  both  above 
and  below  the  current  sheet.  The  observation  point  is  above  the  current  sheet 
when  the  $  angle  is  about  13b0,  ana  below  when  $  -  135°.  The  abrupt  change  in 
angle  signals  the  current  sheet  passage,  or  the  so-called  "sector  boundary 
cross i ng . " 

A.  Outputs 

The  parameters  cnosen  in  Section  2  determine  automatically  the  locations 
of  th«.  earth  an J  a  specific  solar  flare  in  our  coordinate  system,  namely  the 
He. iocentnic  Equatorial  (HF.Q)  coordinate  system,  where  the  zero  longitude  lies 
i'.  ,g  trr  .  ntersection  line  of  the  solar  equatorial  plane  and  the  ecliptic. 
!  •.  th*.  foi  .owing,  we  shall  describe  briefly  each  output.  For  demonstrati  on 
. V„ ;e.  ,  wo  consider  a  hypothetical  flare  which  is  assumed  to  be  reported  by 
pou.  uer  Solar  Observatory,  begining  at  1200  UT  on  December  8,  1990,  the 
I  loo:,  NuO ,  W0G  (namely  at  the  center  of  the  solar  disk',  reaching  the 
m.a  mum  ;  -.tensity  at  1  300  UT ,  and  ending  at  1900  UT:  the  intensity  is  2B. 
0UTPUT_ 3 

Th.s  output  snows  the  projected  location  of  a  hypothetical  solar 
flare  oy  a  circle  on  the  Stanford  source  surface  map  in  Carrington 

1  8 
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long;  tuue- latitude  coordinates  (Output  1 ) .  Since  it  is  assumed  that  a 
flare  takes  place  on  the  center  of  the  solar  disx  in  our  particular 
example,  the  center  of  the  sun,  the  flare  and  the  earth  are  located  on  a 
single  solar  radial  line.  Therefore ,  the  projections  of  both  the  flare 
and  the  earth  on  the  source  surface  coincide  at  1200  UT,  December  8,  1990 
(Output  3).  The  observed  location  of  the  solar  disk  is  automatically 
converted  to  a  point  in  Carrington  coordinates  by  specifying  the  date  ana 
time . 

OUTPUT  aa 

This  output  snows  the  solar  current  sheet  to  a  distance  of  2  AU 
viewed  from  the  longitude  of  the  earth  at  the  time  of  the  solar  flare,  at 
‘200  LIT ,  December  8,  1990.  The  portion  of  the  orbit  of  the  earth  above 
the  sheet  can  be  seen  as  an  elliptical  curve.  From  this  output,  one  can 
see  that  the  earth  is  located  aoove  the  current  at  the  time  of  the 
flare.  The  relative  location  of  the  earth  with  respect  to  the 
heliospheric  current  sheet  is  useful  in  predicting  changes  of  the  IMF 
angle  ij>. 

OUTPUT  Vo 

This  output  is  optional  and  shows  the  solar  current  sheet  to  a 
distance  of  p  AU . 

OUTPUT  ■ 

Equatorial  plane  plots  of  the  disturbed  IMF  are  designed  to  monitor 
the  relative  location  of  the  inter  planetary  disturbance  and  the  earth 
every  six  hours  after  flare  onset  as  the  shock  wave  propagates  outward 
from  the  sun.  The  radius  of  the  plots  extend  to  2  AU.  The  location  of 
the  earth  is  specified  by  the  *  mark.  Space  probe  locations  (if  any)  can 
be  .ndicated  by  PI,  P2,  etc.  If  a  deep  space  probe  at  some  point  between 
trie  sun  and  the  earth  is  available,  the  disturbed  solar  wind  quantities 
it  tne  spacecraft  location  can  be  predicted  and  compared  with  the  space¬ 
craft  observations  to  calibrate  them  well  before  the  expected  disturbance 
arrives  at  the  earth. 

The  f creoaster  will  be  able  to  see  Output  9  on  a  CRT  screen,  so  that 
ne  cm  infer  that  the  shock  wave  will  arrive  at  the  earth  within  a  few 
hour  3  after  the  onset  of  this  particular  flare.  In  this  particular 
exams: ",  the  flare  occurred  at  the  center  of  the  solar  disk  (central 
meridian) ,  so  that  the  center  line  of  the  shock  wave  is  directed  very 
closely  along  the  sun-earth  line;  note  that  the  earth  moved  by  -1.5° 
if  ter  flare  onset.  The  forecaster  can  see  that  the  shock  wave  will 
urr"lv<’  at  the  earth  between  00  and  06  UT  on  December  10  in  this 
nypot  r.'Ot  ica  1  example . 

OUTPUT  u 

A  f . ve-day  forecast  of  sciar  wind  parameters  is  generated  at  the 
earth’s  location  as  a  function  of  time.  The  following  parameters  are 
plotter.;  V  sola.-  wine,  bulk  speed  (km/sec) ;  N  =  solar  wind  number 
tens . ty  ,  arm ) ;  3  -  IMF  magnitude  (n T);  G  =  tne  IMF  latitude  angle 
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Output  5 


(angle  the  IMF  vector  makes  with  the  x-y  plane)  (degrees);  $  =  the 
azimuthal  angle  of  the  IMF  vector  (degrees). 

The  forecaster  will  be  able  to  see  Output  6  on  a  CHT  screen  and  can 
predict  not  only  the  exact  arrival  time,  but  also  how  the  quantities  V, 
N ,  B,  0,  4> )  will  change  as  a  function  of  time.  By  the  time  the  shock 
reaches  the  earth,  it  develops  into  the  forward  and  reverse  shocks  in 
this  particular  case,  as  one  can  see  from  a  double  peak  of  the  speed. 
The  impulsive  increase  of  both  the  density  and  the  IMF  magnitude  B  are 
confined  in  the  two  shocks.  In  this  particular  example,  changes  of  the 
IMF  angle  0  are  negative  and  large.  The  IMF  angle  <p  changes  from  135°  to 
315°,  indicating  the  passage  of  the  heliospheric  currents.  In  general, 
changes  of  the  angle  0  are  much  more  complex  and  depend  greatly  on  the 
relative  location  between  the  flare  and  the  earth  (longitude  and  lati¬ 
tude  ) . 

OUTPUT  6a 

Plots  (V,  N,  B,  Q,  4> )  are  produced  at  specified  space-probe 
locations  between  the  sun  and  the  earth.  This  particular  output  allows 
us  to  calibrate  results  from  the  solar  wind  code  by  comparing  the 
predicted  values  with  observed  values.  For  example,  if  a  Venus  orbiter 
is  located  near  the  sun-earth  live,  this  output  may  be  useful. 

OUTPUT  6b 

Plots  (V,  N,  Bx,  By,  Bz)  are  also  produced  at  specified  space-probe 
Locations  between  the  sun  and  the  earth  in  order  to  calibrate  the  solar 
wind  code  oy  comparing  the  predicted  parameters  with  observed  values. 

OUTPUT  7 

This  output  is  designed  for  the  forecaster  to  predict  how  the 
expected  geomagnetic  storm  develops  as  a  function  of  time.  The  following 
quantities  are  plotted  as  a  function  of  time  at  the  Earth's  location: 

(1)  e  =  the  solar  wind-magnetosphere  dynamo  power. 

(2)  4>p  -  the  total  potential  drop  across  the  polar  cap. 

(3)  Dst  =  the  average  decrease  of  the  horizontal  component  at  the 

equator  (a  measure  of  the  ring  current  intensity). 

(4)  AE  =  the  auroral  electrojet  index. 

(5)  Uj  =  the  total  energy  dissipation  rate  in  the  polar 

ionosphere . 

(6)  ATm  =  an  increased  amount  of  temperature  in  the  auroral 

thermosphere . 

Thr  following  two  quantities  are  also  available  on  demand: 

(7)  Kp  =  the  planetary  (3~hourly)  magnetic  index 

(8)  ap  =  the  linearized  planetary  magnetic  index  T,,. 

Output  7  shows  an  example  of  the  predicted  e,  the  potential  drop 
0pc,  Dst  (CALDST )  ,  AE(CALAE),  Uj-  and  AT*,.  In  this  particular  case,  c 
increases  sharply  during  the  first  six  hours  of  December  10,  1990.  Then 


it  decreases  rather  rapidly  from  -07  to  09  UT  and  slowly  afterward.  The 
storm  'ends'  by  1200  UT  on  December  12.  As  expected,  the  other  quanti¬ 
ties  varied  in  similar  ways. 

Thus,  the  Output  7  allows  the  forecaster  to  predict  how  the  storm 
will  develop  in  terms  of  the  geomagnetic  indices  AE  and  Dst  as  a  function 
of  time.  The  last  two  quantities  are  important  also  in  accessing  an 
increase  of  atmospheric  drag  which  will  be  experienced  by  satellites. 
The  above  quantities  depend  greatly  on  the  IMF  angle  0  which  often  has 
most  complex  changes,  so  that  the  resulting  storm  development  is  most 
often  much  more  complex  than  the  example  shown  here. 

OUTPUT  8 

Output  8  provides  the  geometry  of  the  open  field  line  region  every 
three  hours  after  the  arrival  of  the  shock  wave  at  the  earth  (namely,  the 
storm  sudden  commencement).  If  an  intense  solar  electron  event  is 
expected,  this  output  will  be  useful  in  predicting  the  geometry  of  the 
bombardment  area  by  energetic  solar  electrons.  At  00  and  03  UT  on 
December  10,  the  polar  cap  size  is  normal.  A  slight  asymmetry  with 
respect  to  the  noon-midnight  meridian  is  produced  Dy  the  IMF  B 
component.  At  06  UT,  the  polar  cap  size  increases,  and  the  asymmetry 
becomes  opposite  to  that  at  00  and  03  UT ;  the  IMF  By  component  changes 
the  sign.  The  geometry  of  the  solar  proton  entry  area  is  known  to  be 
appreciably  longer  than  that  of  the  electrons.  Output  7  information 
should  be  available  to  astronauts  in  a  polar  orbiting  shuttle.  It  may  oe 
convenient  to  note  that  the  entry  area  of  solar  protons  is  appreciably 
larger  than  that  of  solar  electrons. 

OUTPUT  9 

Output  9  provides  roughly  the  geometry  of  the  auroral  oval  every 
three  hours  after  the  arrival  of  the  shock  wave  at  the  earth  (namely,  the 
sudden  commencement).  The  size  of  the  oval  is  quantified  in  terms  of  the 
Q  index  (a  quarter-hourly  Kp  index).  Output  9  shows  a  sample  of  the  CRT 
display  oval  during  the  storm.  One  can  see  that  the  oval  is  enlarged 
considerably  at  04  and  07  UT,  but  is  contracted  by  10  UT.  In  the  future, 
the  auroral  distribution  should  be  adjusted  for  the  Dst  index. 


b .  Second  Gener at ion  Schem e 

When  a  large  number  of  codes  are  Integrated  to  develop  a  complex  scheme, 
the  accuracy  and  efficiency  of  the  entire  scheme  is  determined  by  the  weakest 
code  component.  We  list  some  of  the  most  important  improvements  to  be  made  on 
the  basis  ef  future  studies. 
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(A)  IMF  component 


The  weakest  part  in  the  first  generation  scheme  lies  in  properly- 
modeling  variations  of  the  north-south  component  (Bz)  of  the  IMF.  This 
particular  subject  has  fallen  into  the  gap  between  two  disciplines,  solar 
physios  and  magnetospneric  physics,  so  that  there  so  far  have  been  only  a 
few  observational  and  theoretical  papers  to  guide  us  in  developing  our 
code  for  this  particular  purpose. 

Basic-illy,  a  solar  flare  causes  changes  of  the  location  of  the  foot 

point  of  magnetic  field  lines  on  the  source  surface.  Such  changes  will 

be  transmitted  to  i nter planetar y  space  and  cause  changes  of  the  IMF  Bz 
component.  At  present,  our  code,  based  on  this  principle,  adopts  one  of 
the  simplest  changes  of  the  foot  point  of  the  field  lines.  The  so-called 
’magnetic  ilouds'  could  be  another  cause  of  the  IMF  Bz  component. 

Altogether  much  conceptual  development  is  still  needed  in  understanding 
IMF  K.  changes,  as  well  as  associated  code  developments. 

j'n  accomplishing  the  task  described  above,  we  have  to  make  the 
following  efforts  for  individual  solar  events. 

(i)  Examine  solar  observations  in  an  attempt  to  find  quantities 

which  are  crucial  in  influencing  the  IMF  Bz  component. 

(ii)  Model  magnetic  field  variations  in  the  chromosphere  and  corona 
during  solar  flares  in  order  to  infer  the  corresponding 
magnetic  field  variations  on  the  source  surface. 

(lii)  Model  the  propagation  of  the  disturbances  in  (ii)  to  the 
location  of  the  earth.  The  opposite  approach  to  (i)— (ii)-(iii) 
may  also  be  useful.  That  is  to  say,  on  the  basis  of  the 
observed  IMF  Bz  valuations  at  1  AU,  infer  magnetic  field 
variations  on  the  source  surface. 

> B)  Flare  parameters 


Needless  to  say  that  the  present  determi  nation  of  the  flare  para¬ 
meters  i s  a  very  tentative  one.  Much  effort  is  needed  to  find  a  suitable 
/.ilu-'  for  each  parameter.  In  this  respect,  it  is  of  great  interest  to 


note  that  Sheeley  et  al .  (1983)  have  shown  that  coronal  mass  ejections 
(CMEs )  are  very  highly  correlated  with  long  duration  soft  x-ray  events. 
CMEs  are  also  well  correlated  with  interplanetary  shook  waves.  There¬ 
fore,  it  may  well  be  that  ip  is  closely  related  to  the  duration  xD  of 
soft  x-ray  events  (say,  ip  -2xp);  both  Xp  and  the  duration  of  soft  x-ray 
events  are  much  longer  than  the  duration  of  optical  flares.  The  so- 
called  "halo-type"  CMEs  observed  by  the  SOLWIND  need  special  attention, 
because  they  indicate  that  a  large-scale  interpl anetary  disturbance  is 
directed  toward  the  the  earth  (Howard  et  al.,  1982).  An  empirical 
formula  for  timing  the  arrival  of  a  solar  flare-initiated  shock  wave  has 
recently  been  determined  by  Smart  and  Shea  (1985).  Their  results  should 
also  be  incorporated  in  our  determination  of  Vp  in  the  future. 

(C)  Magnetosphere- ionosphere  coupling 

The  second  weakest  part  of  our  scheme  is  that  we  rely  on  several 
empirical  relationships  in  the  second  part  of  the  code.  They  are: 

(i)  Determination  of  the  polar  cap  potential  <J>  by  e  or  AE. 

PL 

(ii)  Determination  of  the  AE  index  by  e. 

These  relationships  should  be  replaced  by  theoretical  ones  when  they 
become  available  in  the  near  future. 

( D )  Ring  c urrent  format  ion 

In  this  part  of  the  scheme,  we  also  use  a  semi-empirical  relation¬ 
ship  in  that  it  is  assumed  that  70?  of  the  power  e  is  consumed  from  the 
ring  current.  This  empirical  relationship  should  be  replaced  by  a 
theoretically-based  relationship  which  use  e  (see  Lee  et  al .  (1983)1. 


^  V*r-S  1 


b.  Concluding  Remarks 


It  is  easy  to  dismiss  the  present  prediction  effort  by  stating  that  such 
an  attempt  is  premature  at  the  present  time.  On  the  other  hand,  such  an 
effort  can  focus  our  attention  to  some  of  the  most  important  issues  in  solar — 
terrestrial  physics,  although  they  may  not  necessarily  be  the  most  important 
issues  in  solar  physics  or  magnetospheric  physics.  As  mentioned  earlier,  the 
weakest  part  of  our  prediction  scheme  arises  from  the  fact  that  causes  of  IMF 
changes  are  not  well  understood.  Thus,  these  causes  are  one  of  the  most 
important  unsolved  problems  in  solar-terrestrial  physics .  Thus,  there  is 
absolutely  no  reason  to  leave  the  most  important  unsolved  problem  unattended 
by  saying  that  such  efforts  are  premature. 

The  first  scheme  is  nothing  more  than  an  attempt  to  demonstrate  the  basic 
principle  of  the  modern  prediction  scheme.  We  would  be  more  than  happy  that 
this  report  could  serve  in  calling  attention  to  some  of  the  most  crucial 
unsolved  problems  in  solar — terrestrial  physics. 

Acknowledgements:  This  work  has  been  funded  by  Air  Force  Contract  No. 

^ 9' 2S-82-K-0035. 


32 


References 


Ahn,  B.-H.,  S.-I.  Akasofu,  Y.  Kamide,  and  J.  H.  King,  Cross-polar  cap 

potential  drop  and  the  energy  coupling  function,  J.  Geophys.  Res.,  89 , 
H028,  1984. 

Akasofu,  S.-I.,  Energy  coupling  between  the  solar  wind  ana  the  magnetos ph -'re, 
Space  Sci .  Rev . ,  28 ,  121,  1981. 

Akasofu,  S.-I.,  Prediction  of  development  of  geomagnetic  storms  by  u.'.ng  t'.- 
solar’  wind-magnetosphere  energy  coupling  function  e,  Planet.  Sp<iC>.. 

29,  'Ibi,  1982. 

Akasofu,  S.-I.,  and  3.  Chapman,  Sol  ar-Ter  res  trial  Physics,  Oxi'ord  Univ.  Pr-’ss, 

1  972 . 

Akasofu,  S.-I.,  and  M .  Roederer,  Dependence  of  the  polar  cap  geometry  or  the 
IMF,  Planet.  Space  Sc: . ,  32,  111,  1984. 

Akasofu,  S.-I.,  K.  Hakamada,  and  C.  Fry,  Solar  wind  disturbances  cvusen  by 
solar  flares:  Equatorial  plane.  Planet .  Space  Sci  . ,  81  ,  ’438,  1  'a 8 3  • 

D’l'ston,  C.,  M.  Dryer,  S.  M.  Han,  and  S.  T.  Wu,  Spatial  structure  of  flare- 
associated  perturbations  in  the  solar  wind  simulated  by  the  two- 
dimensional  numerical  MHD  model,  J ,  Geophys .  Res . ,  86 ,  828,  1  981. 

Fry,  C.  G.,  Three-dimensional  structure  of  the  heliosphere:  yuiet-time  and 
disturbed  periods,  Ph.D.  Thesis,  University  of  Alaska,  Fairbanks,  (344 
pages)  May  1985. 

Gislason,  G. ,  M.  Dryer,  Z.  K.  Smith,  S.  T.  Wu,  and  S.  M.  Hon,  Interplanetary 
disturbances  produced  by  a  simulated  solar  flare  and  equatoriaily- 
fluctuating  heliospheric  current  sheet,  Astrophys.  Space  Sci.,  9 8 ,  '49, 
1984. 


33 


Hakamada,  K .  ,  and  S.-I.  Akasofu,  Simulation  of  three-dimensional  solar  wind 
disturbances  and  resulting  geomagnetic  storms,  Space  Sci.  Rev.,  3J_,  3, 

1  982. 

Howard,  R.  A.,  D.  J.  Michels,  N.  R.  Sheeley,  Jr.,  and  M.  J.  Koomen,  The 
observation  of  a  coronal  transient  directed  at  earth,  Ap ■  J . ,  263.  1  01, 

1  982. 

Lee,  L.  C.,  G.  Corrick,  and  S.-I.  Akasofu,  On  the  ring  current  energy 
injection  rate.  Planet.  Space  Sci.,  31 »  901,  1983. 

Olmsted,  C.  D .  ,  and  S.-I.  Akasofu,  One-dimensional  kinematics  of  particle 
stream  flow  with  application  to  solar  wind  simulations,  Planet .  Space 
Sol .  (in  press),  1985. 

Perreault,  P.,  and  S.-I.  Akasofu,  Study  of  geomagnetic  storms,  Geophys.  J. 
Roy .  Astron.  Soc . ,  54 ,  547,  1978. 

PrGiss,  G.  v»’.,  Correlation  between  upper  atmospheric  temperature  and  solar 
wind  conditions,  J.  Geophys.  Res,  (in  press),  1985. 

Re.*':',  ?.  H.,  h.  W.  Spiro,  and  T.  W.  Hill,  Dependence  of  polar  cap  potential 
drop  on  i nter planetary  parameters,  J,  Geophys.  Res,,  86,  7639,  1981. 

Greeley,  N .  h.,  Jr.,  R.  A.  Howard,  M.  J.  Koomen,  and  D.  J.  Michels, 
Associations  between  coronal  mass  ejections  and  soft  x-ray  events,  Ap . 
J^,  2_72,  3«9,  1983. 

Smart ,  D.  F.,  and  M.  A.  Shea,  A  simplified  model  for  timing  the  arrival  of 
solar  flare-ini tiated  shocks,  J.  Geophys.  Res,,  90,  183,  1985. 


34 


